EPJ manuscript No. 

(will be inserted by the editor) 



Ground-state correlation properties of charged bosons trapped in 
strongly anisotropic harmonic potentials 



o 

(N 
S3 

a 

(N 



Przemyslaw Koscik 1 , Anna Okopinska 1 

Institute of Physics, Jan Kochanowski University 
ul. Swi§tokrzyska 15, 25-406 Kielce, Poland 

the date of receipt and acceptance should be inserted later 

Abstract. We study systems of a few charged bosons contained within a strongly anisotropic harmonic trap. 
A detailed examination of the ground-state correlation properties of two-, three-, and four-particle systems 
is carried out within the framework of the single-mode approximation of the transverse components. The 
linear correlation entropy of the quasi-lD systems is discussed in dependence on the confinement anisotropy 
and compared with a strictly ID limit. Only at weak interaction the correlation properties depend strongly 
on the anisotropy parameter. 
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1 Introduction 



Systems of few interacting particles confined by an exter- 
nal potential are of increasing interest in view of their ap- 
plication to model various nanostructures fabricated in an 
artificial way. The Schrodinger equation for harmonically 
trapped particles which interact through a Coulomb po- 
tential can be used to simulate different physical systems, 
such as semiconductor quantum dots [1] or electromagnct- 
ically trapped ions [2]. Besides of fermionic systems, also 
the bosonic ones have been thus considered in various the- 
oretical contexts both in the 3D [31415] and 2D [5] case. 
The properties of quasi-lD systems are, however, rather 
rarely studied, except for the models of bosons with a con- 
tact potential |7l8l9ll0lllll2ll3j . In this work, we discuss 
the case of N identically charged spinless bosons confined 
in a quasi-lD trap. Such a trapping potential may be ex- 
perimentally realized using highly anisotropic harmonic 
traps where the radial confinement is much tighter than 
the axial one. 

We discuss the effects of both the number of particles 
and the control parameters of the system on the correla- 
tion properties. We consider various characteristics such as 
single-particle reduced density matrix, single-particle den- 
sity, and linear entropy. In particular, the influence of the 
anisotropy on the correlations within a quasi- ID structure 
is investigated. 

Our paper is organized as follows. In Section [5] we 
present the model and provide an analytical formula for 
the effective interaction potential. Section [3] surveys the 
quantities we use to characterize correlations in the sys- 
tem. In section [4] the limit of e — > oo is discussed. The 
results are presented in Section [5] and a summary of our 
conclusions is given in Section [5] 



2 Effective Hamiltonian 

Consider a system of N Coulombically interacting parti- 
cles trapped in a 3D axially-symmetric harmonic potential 
with a Hamiltonian given by 
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where pi = y/y? + zf. After the scaling r i-> J ^£^r, 
E i— > hu x E, the Schrodinger equation takes the form 



ff#(ri, r 2 , r N ) = EW (n, r 2 , r N ), 



with 
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The dimensionless coupling g 



1 



^p- represents the 

ratio of the Coulomb interaction to the longitudinal trap- 
ping energy scale and the dimensionless parameter e = 
measures the anisotropy of the trap. We focus our atten- 
tion on the strong anisotropy case, e> 1, when the sys- 
tem becomes quasi-lD. In this case the particles may be 
assumed to stay in the lowest energy state of the tranverse 

Hamiltonian H± = — + ^e 2 p 2 and the one-mode ap- 
proximation is justified. The N— body wave function may 
be taken in the form 

!?(ri,r 2 ,...,rjv) = i'(xi, x 2 , x N )njl 1 ip(y i )(p(z. l ), (4) 
i „ 2 

where ip(z) = (-^)*e~^~ and tp is assumed to be a real 
function. After substituting (jU) into ©, multiplying it 
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from the left by (p(yi)<p(zi)...<p(yN)ip(zN) and integrating 
over 2/1,2/2, -,2/iV, z 1} z 2 , ...,z N we arrive at 

H 1D ip(x 1 ,X2,—,x N ) = E 1D ip(xi,x 2 ,...,x N ), (5) 

where the quasi-lD Hamiltonian has a form 
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Hid = Y}-^d^ + 2^ ]+ ^ 9UlD{ - Xt,X]) + Ne -^ ] 

i—l 1 i<j 

The effective interaction potential 

[<p(yi)«p{y2)v(zi)ip(z 2 )] 2 



Ui D (x 1 ,x 2 ) 



|ri - r 2 
is calculated to be given by 

/ 67T g(x 2 -x 1 ) 



-dy\dy 2 dzxdz{7) 



U 1D {xi,x 2 ) 



(l-erf[\/-\x 2 - Xl \}), (8) 



where erf(z) is the error function. In Fig. [T] the effective 
potential (J5J is compared with the pure Coulomb potential 
for different values of e. We can notice that the larger is 
the value of e, the closer to the origin does the effective 
potential begin to exhibit Coulomb behaviour. 




Fig. 1. Effective interaction potential (JSJ) as a function 
of x = x 2 — x\ for the anisotropy parameter e = 30, 100 
compared with the Coulomb potential x~ x . 



We have tested the applicability of the single mode ap- 
proximation for the two-particle system (N = 2) by com- 
paring the ground-state energy E\o obtained from the 
ID Hamiltonian ([5]) with the energy E determined from 
the full 3D Hamiltonian ((3]). The relative error AE = 
(\E — Eiu\)/Eq as a function of e is shown in Fig. [2] for 
different values of the coupling constant. One can con- 
clude that anisotropy ratios e > 5 are sufficiently large for 
employing the single mode approximation. 
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Fig. 2. Relative ground-state energy error 

AE = \E - E XD \/E as a function of e for N = 2 
and g = 1,5, 20. 



for spinless particles defined as 

p(r,r')= ^(r,r 2 ,...,r A r)if(r',r 2 ,...,r Ar )d 3 r2...d 3 r w . 

(9) 

In the one-mode approximation (0| the RDM takes the 
form 

p(r,r ) = <p(y)<p(y )<p(z)(p(z )p 1D {x,x ), (10) 
with the ID effective RDM given by 



Pid(x,x ) = I ... I ip(x,x 2 ,...,x N )ip(x ,x 2 ,...,x N )dx 2 ...dx N . 

(11) 

The effective RDM can be represented in the Schmidt form 

DO 

pi D (x,x ) = y^ j \ivi(x)vi{x ), (12) 
;=o 

where {vi(x)} are the natural orbitals and their occupan- 
cies {A;}. We will concentrate on discussing the linear en- 
tropy 

r r 

Pid(x, x') 2 dx'dx. 



L — I — 



(13) 



3 Correlation characteristics 

A basic tool to investigate two-body correlations in the 
system is the one-particle reduced density matrix (RDM) [2] 



which can be expressed in terms of A; as L = 1 — X 2 . 
It gives indication of the spread of terms in the Schmidt 
decomposition (|T2l and is one of the popular measures of 
correlation |15ll6ll7ll8j . 



4 Strictly ID limit 

In the strictly ID limit of e — > oo, the interaction potential 
of the considered system, g/\xj — Xj\, diverges at Xi = Xj 
for any finite g. This causes divergences at short distances 
in calculating the energy of bosonic systems. Usually those 
divergences are cured by performing calculation for an 
anisotropic 3D system with finite e and observing the lim- 
iting behavior at e — > oo. However, as noticed recently |19) . 
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the calculation may be performed directly for the ID sys- 
tem since its ground-state wavefunction can be related via 
Bose-Fermi mapping to the lowest energy antisymmetric 9 ' 
TV-particle wavefunction ipp as 



ip(x 1 ,X2,..,x N ) = \ip F (xi,x 2 ,..,x N )\. 



(14) 



Therefore, the system of bosons gets fermionized and the 
ultraviolet divergencies are cured in a natural way when 
determining ipp by the standard configuration interaction 
method. 



5 Results and discussions 

Wc shall consider the ground-state of many-body systems 
with the number of particles TV = 2, 3 and 4. In order 
to reveal qualitatively the nature of correlations in the 
quasi- ID limit, we first discuss the case of large anisotropy, 
e = 30. The ground-state TV-particle wave function of the 
Hamiltonian ([SJ is calculated with the quantum diffusion 
algorithm [20] and used to determine p(x, x ) by numerical 
integration of Eq. (fTTj) . 
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Fig. 3. Gray-scale plots of the effective RDM for systems 
of TV = 2,3 and 4 particles with confinement anisotropy 
e = 30 at various interaction strengths g. 



Grey-scale plots of p(x, x ) for the considered systems 
are shown in Fig. [3] at three values of the interaction 
strength g. As one can see, the off-diagonal elements of the 
RDM diminish with increasing g, which indicates a loss of 
spatial coherence in the system. This is due to the re- 
pulsive nature of the interaction encouraging localization 
of the particles, which tend to separate from each other. 
A clear deviation from the circular structure p(x, x ) = 




TV = 3 



TV = 4 



♦ 
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Fig. 4. Gray-scale plots of the effective RDM for systems 
of TV = 2,3 and 4 particles with confinement anisotropy 
e = oo at various interaction strengths g. 



of the noninteracting case (g = 0) is clearly 
1. Interestingly enough, the re- 



7r 2 e 

observed already at g 
suits of Fig. [3] show that the onset of crystallization shows 
up at g = 5, independently on the number of particles. 




5 = 2 





Fig. 5. Single-particle density of TV = 3 particles confined 
with anisotropy e = 30 (dashed curve) and e = oo for 
various interaction strengths g. 



The results at finite anisotropy e = 30 may be com- 
pared with those obtained in the strictly ID limit that are 
given in Fig. [4] In this case we determined tpp, and thereby 
tpB, using the standard configuration interaction method 
based on harmonically trapped single-particle cigenfunc- 
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tions (fin . Comparing the results of Fig. [3] with the ones 
of Fig. |4j one can notice that the anisotropy parameter 
influences the behaviour of the RDM only in the regime 
of small values of g. As a matter of fact, already at g = 5 
the RDM calculated at e = 30 reproduces quite well the 
one calculated in the strictly ID limit, regardless of the 
number of particles. We observe that below this value an 
increase in the anisotropy parameter has the effect of an 
increase in the interparticles distances. Otherwise stated, 
the onset of the Wigner crystalization in the strictly ID 
limit appears at value smaller than that in the case of 
e = 30. To make the above more clear, we present in Fig. 
[5] the single particle densities of the 3-particle system, cor- 
responding to the N — 3 results of Figs. [3] and 2) 

Next, we explore the effect of the anisotropy parameter 
e on the correlation properties of systems containing N = 
2, 3 and 4 bosonic particles. In Fig. [5] the ground-state 
linear entropy L at different values of e is compared with 
the result obtained in the fermionizcd strictly ID limit of 
e — > oo. The results are presented as a function of the 
dimensionless parameter g. 

We see that in each considered case, the linear entropy 
increases with g and above g f=s 5 saturates at a constant 
value that is insensitive to the anisotropy parameter e. 
This may be attributed to the fact that at large g the av- 
erage distances between the particles are large enough so 
that the effective interaction exhibits pure Coulomb be- 
haviour. The values of g at which saturation takes place 
seem to increase slightly with N. A look at Fig(3] allows 
us to conclude that they coincide roughly with the critical 
values of crystallization. As can be inferred from Fig. [5J 
the linear entropy L increases with the number of parti- 
cles in the system, with the increase getting smaller for 
larger N. The effect becomes less pronounced at weaker 
interactions and disappears in the limit of g — > 0, when 
L —> 0, regardless of the number of particles. At small val- 
ues of g, the entropy L depends strongly on the anisotropy, 
being larger at larger e. Again, this can be qualitatively 
understood by referring to the distances of the particles, 
namely they are small at small values of g and are thus 
in the regime where the effective interaction potential (O 
strongly depends on e. 

The entropy in the strictly ID limit of e — > oo has been 
calculated by taking into account the Bose- Fermi mapping 
(fl"4")) . In the limit of g 0, the wavefunction ipp reduces 
to a single Slater determinant and (fT4)l takes the form 
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ip(x 1 ,x 2 , .-,x N ) 



1 



\det^L N =M°{x 3 ))\- (15) 
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Fig. 6. Linear entropy in the ground-state of the system 
of © for N — 2, 3 and 4 particles as a function of g. 
Full curve, e = 5; dashed curve, e = 30; dot-dashed curve, 
e = 100. The red line is the linear entropy in the strictly 
ID limit. 



approximation which becomes exact in this limit |18j . The 
calculation analogous to that performed by us in the 2D 
case [T5] gives 



The values of Lirj are calculated to be about 0.36, 0.51 
and 0.6 for TV = 2, 3 and 4, respectively. The occurrence 
of fermionization for g ^ results in the discontinuity of 
the linear entropy in the point g = 0. Importantly, we 
observe that the smaller is g and/or N the larger is the 
anisotropy parameter at which the linear entropy of the 
quasi- ID system reaches the fermionic behaviour of the 
strictly ID gas (e — > oo). 

The asymptotic value of Liu at g — > oo can be calcu- 
lated analytically in the case of A = 2, using the harmonic 



r9-s-oo 



= 1 - 



-- + V3w 0.518, (16) 



which is in agreement with our numerical result. For N = 
3 and N = 4 the values of Lf~^°° are found numerically 
to be about 0.68 and 0.77, respectively. To fully reveal the 
correlation effects in the strong interaction limit, it would 
be desirable to obtain the full dependence of L 9 ^ 00 on N. 
This is a topic for future investigation. 
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6 Summary 

We investigated the ground-state correlation properties of 
the systems composed of two, three, and four particles 
in a strongly anisotropic harmonic trap. Within the one- 
mode approximation, we studied the influence of both the 
number of particles N and the anisotropy parameter e 
on the correlation properties of the systems in the whole 
range of repulsive interaction strength g. As a general 
trend we found that the linear entropy L increases with 
increasing g. At small g (weak interaction and/or strong 
confinement), the entropy of the considered systems de- 
pends heavily on the anisotropy of the trap, being larger 
at higher anisotropy. Linear entropy is the largest in the 
limit of e — > oo, when fermionization takes place for any 
g =/= 0. At large g (strong interaction and/or weak confine- 
ment), the entropy L saturates at a value that does not 
depend on e and is greater the larger is N. The value of 
g at which saturation takes place hardly depends on the 
number of particles, shifting only slightly towards larger 
values with increasing N. 

The practical realizations of the model discussed in our 
work can be achieved experimentally in linear ion traps 
where the confining forces in the longitidunal direction 
are much softer than the radial confinement. In th e case 
of singly charged ions the parameter g — ke 2 ^ J^z , where 

k is the Coulomb constant and m is the ion's mass, can 
be controlled by the axial trapping frequency u> x . In cur- 
rent experiments the values of oj x are below 1MHz, which 
corresponds to the crystalline phase with micrometer dis- 
tances between the ions [21] and a value of g larger than 
10 6 . In this regime, the correlation depends very weakly 
on anisotropy and the strictly one-dimensional approxi- 
mation works well even at moderate values of e. 

We think a more extensive analysis of the effects of 
anisotropy for a larger number of particles is desirable to 
get a deeper insight into the properties of strongly corre- 
lated bosons with Coulomb interaction. 
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